Optimally solving a transportation problem using 
Voronoi diagrams^ 



Darius GeiB a , Rolf Klein a , Rainer Penninger a '*, Giinter Rote b 

Rheinische Friedrich- Wilhelms Universitat Bonn, Institute of Computer Science I, 

Friedrich-Ebert-Allee 144, D-53113 Bonn, Germany. 
b Freie Universitat Berlin, Institut fur Informatik, Takustrafle 9, D- 14195 Berlin, 

Germany. 



Abstract 

We consider the following variant of the well-known Monge-Kantorovich 
transportation problem. Let S be a set of n point sites in IR d . A bounded 
set C C M. d is to be distributed among the sites p G S such that (i), each 
p receives a subset C p of prescribed volume and (ii), the average distance of 
all points z of C from their respective sites p is minimized. In our model, 
volume is quantified by a measure //, and the distance between a site p and 
a point z is given by a function d p (z). Under quite liberal technical assump- 
tions on C and on the functions d p (-) we show that a solution of minimum 
total cost can be obtained by intersecting with C the Voronoi diagram of 
the sites in S, based on the functions d p (-) equipped with suitable additive 
weights. Moreover, this optimum partition is unique, up to sets of measure 
zero. Unlike the deep analytic methods of classical transportation theory, 
our proof is based directly on geometric arguments. 
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1. Introduction 

In 1781, Gaspard Monge [8] raised the following problem. Given two sets 
C and S of equal mass in K d , transport each mass unit of C to a mass unit of 
S at minimal cost. More formally, given two measures /i and v, find a map 
/ satisfying = v(-) that minimizes 



where d(z,z') describes the cost of moving z to z'. 

Because of its obvious relevance to economics, and perhaps due to its 
mathematical challenge, this problem has received a lot of attention. Even 
with the Euclidean distance as cost function d it is not at all clear in which 
cases an optimal map / exists. Progress by Appell [1] was honored with a 
prize by the Academy of Paris in 1887. Kantorovich [7] achieved a break- 
through by solving a relaxed version of Monge's original problem. In 1975, 
he received a Nobel prize in Economics; see Gangbo and McCann [3] for 
mathematical and historical details. 

While usually known as the Monge-Kantorovich transportation problem, 
the minimum cost of a transportation is sometimes called Wasserstein metric 
or, in computer science, earth mover's distance between the two measures /i 
and v. It can be used to measure similarity in image retrieval; see Rubner et 
al. [TO]. If both measures \x and v have finite support, Monge's problem be- 
comes the minimum weight matching problem for complete bipartite graphs, 
where edge weights represent transportation cost; see Rote [Sj, Vaidya [12] 
and Sharathkumar et al. [TT] . 

We are interested in the case where only measure v has finite support. 
More precisely, we assume that a set S of n point sites pi is given, and 
numbers \ > whose sum equals 1. A body C of volume 1 must be split 
into subsets of volume \ in such a way that the total cost of transporting, 
for each i, all points of Cj to their site Pi becomes a minimum. In this setting, 
volume is measured by some measure fi, and transport cost d(z, z') by some 
measure of distance. 

Gangbo and McCann [3] report on the cases where either d(z, z') = 
h(z — z') with a strictly convex function h, or d(z,z') = l(\z — z'\) with 
a non-negative, strictly concave function I of the Euclidean norm. As a con- 
sequence of deep results on the general Monge-Kantorovich problem, they 
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prove a surprising fact. The minimum cost partition of C is given by the 
additively weighted Voronoi diagram of the sites pi, based on cost function d 
and additive weights u>j. In this structure, the Voronoi region of pi contains 
all points z satisfying 

d(pi, z) — Wi < d(pj, z) — Wj for all j ^ i. 

Villani [13] has more recently observed that this holds even for more general 
cost functions, that need not be invariant under translations, and in the case 
where both distributions are continuous. 

Figure [T] depicts how to obtain this structure in dimension d = 2. For 
each point site p G S we construct in IR 3 the cone {(z, d(p, z) —w p ) \ z G 1R 2 }. 
The lower envelope of the cones, projected onto the W-plane, results in 
the Voronoi diagram. Independently, Aurenhammer et al. [2] have studied 




Figure 1: An additively weighted Voronoi diagram as the lower envelope of cones. 



the case where d(z,z') = \z — z'\ 2 . Here, a power diagram (see [3]) gives 
an optimum splitting of C. In addition to structural results, they provide 
algorithms for computing the weights Wj. Regarding the case where the 
transportation cost from point z to site pi is given by individual cost functions 
d Pi (z), it is still unknown how to compute the weights W{. 

In this paper we consider the situation where the cost d(pi, z) of trans- 
porting point z to site Pi is given by individual distance functions d Pi (z). We 
require that the weighted Voronoi diagram based on the functions d Pi (-) is 
well-behaved, in the following sense. Bisectors are (d — l)-dimensional, and 
increasing weight Wi causes the bisectors of pi to sweep (i-space in a contin- 
uous way. These requirements are fulfilled for the Euclidean metric and, at 
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least in dimension 2, if the distance function d p (z) = d(z — p) assigned to 
each site p is a translate of the same strictly convex distance function d(-). 

We show that these assumptions are strong enough to obtain the result 
of [U 12]: The weighted Voronoi diagram based on the functions d p .(-) op- 
timally solves the transportation problem, if weights are suitably chosen. 
Whereas [3] derives this fact from a more general measure-theoretic theo- 
rem, our proof uses the minimization of a quadratic objective function and 
geometric arguments. The purpose of our paper is to show that such a simple 
proof is possible. 

After stating some definitions in Section[2]we generalize arguments from [2] 
to prove, in Section [3j that C can be split into parts of arbitrary given vol- 
umes by a weighted Voronoi diagram, for a suitable choice of weights. Then, 
in Section |4j we show that such a partition is optimal, and unique. In Sec- 
tion [5] we show that if C is connected, these weights are uniquely determined, 
up to addition of a constant. 

2. Definitions 

Let p be a measure defined for all Lebesgue-measurable subsets of M. d . 
We assume that p and the (i-dimensional Lebesgue measure are mutually 
continuous, that is, p vanishes exactly on the sets of Lebesgue measure zero. 

Let S denote a set of n point sites in M d . For each p G S we are given a 
continuous function 

d„ : — y K>o 

that assigns to each point z of M. d a nonnegative value d p (z) as the "distance" 
from site p to z. 

For p 7^ q G S and 7 G R we define the bisector 

B 7 (p, q) := { z G R d I dp{z) - d q {z) = 7 } 

and the region 

R 7 (p } q ) : = { z e W l I d p (z) - d q (z) < 7 }. 

The sets i? 7 (p, q) are open and increase with 7. Now let us assume that for 
each site p G S an additive weight w p is given, and let 

W = (Wx,W2, ...,W n ) 
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denote the vector of all weights, according to some ordering pi,p2, ... of 5". 
Then, B Wi _ Wj (pi,pj) is called the additively weighted bisector of pi and pj, 
and 

VR w (pi,S) := f]R w% - Wj (Pi,Pj) 

is the additively weighted Voronoi region of Pi with respect to S. It consists 
of all points z for which d Pi (z) — wi is smaller than all values d P] (z) — Wj, by 
definition. As usual, 

V w (S):=R d \{JVR w (p t ,S) 

i 

is called the additively weighted Voronoi diagram of S; see [3J. Clearly, 
VR w (Pi, S) and V W (S) do not change if the same constant is added to all 
weights in w. Increasing a single value Wi will increase the size of p^s Voronoi 
region. 




(i) («) 

Figure 2: (i) Sets i? 7 (p, q) for increasing values of 7. (ii) An additively weighted Voronoi 
diagram. 

Example. Let d = 2 and <i p (z) = \p — z\ the Euclidean distance. Given 
weights w p ,Wg, the bisector B Wp - Wq (p,q) is a line if w p = w q , and a branch 
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of a hyperbola otherwise. Figure j2](i) shows how B 1 (p,q) sweeps across the 
plane as 7 grows from —00 to 00. In fact, when |7| = \p — q\, the bisector 
degenerates to a ray, and for larger |7|, the bisectors are empty. The bisector 
B y (p, q) forms the boundary of the set R~,(p, q), which increases with 7. Each 
bounded set C in the plane will be reached at some point. From then on 
the volume of C nR 7 (p, q) is continuously growing until all of C is contained 
in i? 7 (p, q). Given n points Pj with additive weights Wj, raising the value of 
a single weight Wi will cause all sets R Wi - Wj (pi,pj) to grow until C is fully 
contained in the Voronoi region V w (pi, S) of pi. 

Figure [2](ii) shows an additively weighted Voronoi diagram V W (S) based 
on the Euclidean distance. It partitions the plane into 5 two-dimensional 
cells (Voronoi regions) , and consists of 9 cells of dimension 1 (Voronoi edges) 
and of 5 cells of dimension (Voronoi vertices). Each cell is homeomorphic 
to an open sphere of appropriate dimension. 

The next definition generalizes the above properties to the setting used 
in this paper. 

Definition 1. A system of continuous distance functions d p (-), where p G S, 
is called admissible if for all p ^ q G S , and for each bounded open set 
C C WL d , there exist values m pq and M pq such that 7 h-> fi (C C\R 7 (p,q)) is 
continuously increasing from to p,{C) as 7 grows from m pq to M pq . Fur- 
thermore, C fl R~,(p,q) = ifj< rn pq and C C R~,(p, q) if M pq < 7. 

By symmetry we can assume w. 1. o. g. m qp = —M pq and M qp = —m pq \ see 
Figure |2](i). We will need the following structural property, which essentially 
says that bisectors have measure 0. 

Lemma 1. For a system of admissible distance functions d p (-), where p G S , 
for any two points p 7^ q G S and any 7 G and for any bounded open set 
C C W l , we have fj, (C D B 7 (p, q)) = 0. 

Proof. By definition of R^(p,q) and B y (p,q), for all e > and 7 G M. the 
following inequality holds: 

n(C n R j+£ (p, q)) > n{C n R^ip, q)) + n{C n B 7 (p, q )) 

We let e decrease to 0. Since, according to Definition [TJ the function 7 1— > 
fi (C fl Ry(p, q)) is continuous, we get 

lim (ji(C n Ry +£ (p, q))) = H{C n Ry(p, q)). 
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This implies /x(C D B y (p, q)) = 0. 



□ 



3. Partitions of prescribed size 

The following theorem shows that we can use an additively weighted 
Voronoi diagram based on an admissible system of distance functions d p to 
partition C into subsets of prescribed sizes. 

Theorem 1. Let n > 2 and let d Pi (-), 1 < i < n, be an admissible system 
as in Definition^ Let C be a bounded open subset o/R d . Suppose we are 
given n real numbers A, > with Ai + A2 + • • • + A n = //(C). Then there is a 
weight vector w = (wi, w 2 , ■ ■ ■ , w n ) such that 

//(C n VR w (p i , S)) = A, 

holds for 1 < i < n. 

If, moreover, C is pathwise connected then w is unique up to addition of 
a constant to all Wi, and the parts Ci = C D VK w (pi, S) of this partition are 
unique. 

Proof. The function 

n 

$H := ]T(MC n VR w (p t , s)) - a,) 2 

measures how far w is from fulfilling the theorem. Since each function 7 y 
fj,(C n Ry(pi,pj)) is continuous by Definition [TJ and because 

\n(c n VR^fe, S)) - pl{c n VR^fo, < 

Wi — Wj 

we conclude that the function $ is continuous, too. 

Let D := max{ M pg \ P 7^ q} with M pg as in Definition [lj Note that 
m pg = —M qp and m P5 < M pq together imply that D > is a positive constant. 
On the compact set [0, D] n , the function $ attains its minimum value at some 
argument w. If $(to) = we are done. Suppose that > 0. Since the 

volumes of the Voronoi regions inside C add up to //(C), there must be some 
sites Pj whose Voronoi regions have too large an intersection with C i.e., of 
volume > Xj, while other region's intersections with C are too small. 
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We now show that by decreasing the weights of some points of S we can 
decrease the value of $. For any point pi G S we consider the continuous 
"excess" function 

4> Pi (w) := /x(c n VR w (Pi, s)) - \i 

that indicates by which amount of volume the region of pt is too large or too 
small. We have 

n 

£> ft («0 = o. 

i=i 

Let T C S be the set of points pi for which <f> Pi (w) attains its maximum 
value, t. By assumption, T is neither empty nor equal to S. We will reduce 
the weights of the points in T by the same amount 5, obtaining a new weight 
vector w'. Note that, by construction, VR w (t, S) D VR w i(t,S) for any site 
t G T, and VR™(s, S) C VR u ,/(s, S 1 ), for any site s E S \ T. 

Let denote the weight vector (u^, . . . ,w' n ) where w[ := w« — 5 if 

Pi G T, and -u^ := otherwise. The volume of the set of points of C, that 
is assigned by V W (S) to some point in T and by V w >(g)(S) to some point in 
S\T is 

Loss(5) :=Y,(Mu)-Mtf(S))) 
teT 

Let t := maxpgs <j> p (w) > and r' := max sg 5\T 4> s {w) < r. Our goal is to 
choose S in such a way that 

Loss(5) = T -^- (1) 

and to show that this choice decreases $ strictly. Since Loss(O) = and 
Loss (5) is a monotone and continuous function, our aim is to solve ([TJ by 
the intermediate value theorem. For all 5 > 2D we observe that Loss(<5) = 
S Pi eT( r +^i) > T holds, since in this case the Voronoi region VR w '(8)(Pi, S) of 
every point pi G T contains no point of C. Now < < r implies that we 
can indeed solve (JTJ by the intermediate value theorem. The latter inequality 
is evident if r' > 0. Otherwise, it follows from — r' < — J2 s es\T 4>s(w) = 
52 seT (j)s(w) = \T\ ■ t < (n - l)r. 

In the following let w' := w'(S), for short. In order to prove the theorem 
we use the following reformulation of the statement <$>(w) > 

Y,(<i>tM 2 - Mw') 2 ) > E i^ w ') 2 - ^H 2 ) ( 2 ) 

teT seS\T 
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2n ' 



Now, in order to prove inequality we give an upper bound on the right- 
hand side of the inequality, which will be compared with a lower bound on the 
left-hand side of the inequality. Let e s > denote the volume increase, in C, 
of the region of site s G S\T. By the choice of 5, we have ^2 seS \ T £ s 
and 

{Mw'f-Mw) 2 ) = ((0 s H + ^) 2 -0,H 2 ) 

s£S\T s£S\T 

ses\T ses\T 
<2r' £ ^ 

s&S\T s£S\T 

I / i\ 2 

,t — r it — t 
2r'— + 



2n V 2n 



We now give a lower bound on the left-hand side of ([2]). If St > denotes 
the volume decrease, in C, of the region of site t G T, then the corresponding 
calculation for the points in T yields: 



2 

T — T / T — T 

2r- 



2n V 2ri 

where we have used the same arguments as above and also that <fit{w) = r, 
by definition, for all teT. It is easy to verify the inequality 

T — T I T — T \ ,T — T I T — T 

It— — > 2r'- 



2n \ 2n I 2n \ 2n 



for t — t' > and n > 2. Hence inequality (|2]) holds. 

We have now shown the existence of a weight vector w' with < 
$(u>). If w' belongs to [0, D] n , we are done. Otherwise, suppose the maximum 
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weight w' m occurs for some point p' m . We add D — w' m to all weights. This 
does not change the Voronoi diagram defined by w', but now the maximum 
weight w' m is exactly D. Next we observe that if now the weight of some point 
Pj is < 0, then its Voronoi region is empty, because C is contained in the 
set R w i -w'XpmiPj), by definition of D. Assigning weight to each of those 
points also results in empty Voronoi regions for those points. Altogether, we 
now have < w\ < D for all pi G S. 

We have shown the existence of a weight vector w' G [0, D] n that satisfies 
$(«/) < $(u>). This contradicts the minimality of Hence the minimum 

value of $ is 0. 

Uniqueness of the partition will be discussed in Section [5j □ 
4. Optimality 

Again, let d p (-), p G S = {pi,P2, ■ ■ ■ ,Pn}, be an admissible system as in 
Definition [TJ and let C denote a bounded and open subset of M. d . Moreover, 
we are given real numbers Aj > whose sum equals n(C). 

By Theorem [I] there exists a weight vector w = (u>i, W2, ■ ■ ■ , w n ) satisfying 

\i{C n VR w (pi, S)) = Xi for i = 1, . . . , n. 

Now we prove that this subdivision of C minimizes the transportation cost, 
i. e., the average distance from each point to its site. It is convenient to 
describe partitions of C by /i-measurable maps /: C — > S where, for each 
p G S, denotes the region assigned to p. Let Fa denote the set of 

those maps / satisfying = Aj for i — 1, . . . , n. 

Theorem 2. The partition of C into regions C{ := Cr\VK w (pi, S) minimizes 

cost(/) := / d f ( z )(z) dfi(z) 
Jc 

over all maps f G F\. Any other partition of minimal cost differs at most 
by sets of measure zero from (Ci)j. 

If the measure \i were a discrete measure concentrated on a finite set 
C of points, the optimum partition problem would turn into the classical 
transportation problem, a special type of network flow problem on a bipartite 
graph. The weights Wi would be obtained as dual variables. The following 
proof mimics the optimality proof for the discrete case. 
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Proof. Let f w be defined by f w {Ci) = {p^} for all i, and / G F A be another 
map. If we define c(z,p) := d p (z), then the cost of map / is 

cost(/) = / (c(zj(z)) -w f(z) +w f{z) ) dfx(z) 
Jc 

(c(z,f(z)) -w m ) dfi(z) + / w m dn(z) (3) 
c Jc 

> / (c(z,f w (z)) -w Mz) ) dfx(z) + / Wf( g )dn(z) (4) 
Jc Jc 

= cost(/ tt )- / w /iu(2) d/i(z) + / w /(2) d/i(z) (5) 

where the inequality between lines ^ and Q follows from the fact that 

c(z, f w (z)) - w /m(2) < c(z, f(z)) - w f ( z ), (6) 

by the definition of additively weighted Voronoi diagrams. Here we have 
assumed for simplicity that the map f w assigns every point p £ V W (S) to 
some weighted nearest neighbor of p in S, according to some tie-break rule. 
Since Lemma [I] implies /x(C D V W (S)) = 0, changing the assignment of those 
points has no influence on cost (/„,). 

Both maps / and f w partition C into regions of volume Aj, 1 < i < n, 
i. e., nif^ipi)) = Kf~ l {Pi)) = for a11 Pi e S- Therefore, 

/ w f{z) dfi(z) = y2 / w f(z) dfi(z) = y2\iWi, 

and the same value is obtained for f w : 

/ WMz) &K Z ) = J2 W M Z ) d ^ z "> = XiWi - 

Now the optimality claim follows from (J5]). 

We still have to prove uniqueness of the solution, up to a set of mea- 
sure 0. Assume that / differs from f w on some set G of positive measure. 
We are done if we show that the inequality between ^ and Q holds as a 
strict inequality. We can remove the points of V W (S) from G, since they 
have measure 0. Then ([6| holds as a strict inequality on G. The dif- 
ference between ^ and Q is then the integral of the positive function 
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g(z) := c(z,f w (z)) - w Mz ) - (c(z,f(z)) - w f{z) ) over a set G of positive 
measure. Such an integral J zeG g(z) dfi(z) has always a positive value: the 
countably many sets { z G G | | > 0(2) > ^-j- } and { z G G \ g(z) > 1 } 
partition G, and thus their measures add up to fi(G). Therefore, at least one 
of these sets has positive measure, and it follows directly that the integral is 
positive. □ 

5. Uniqueness 

In contradistinction to Theorem |2j the weight vector w and the resulting 
weighted Voronoi partition in Theorem [T] is not unique: if C is disconnected, 
changing w may move the bisectors between different connected components 
of C without affecting the partition, or only changing it in boundary points 
of r. 

However, uniqueness can be obtained under the additional assumption 
that C is pathwise connected. Then in Theorem [TJ the weight vector w = 
(wx, ■ ■ ■ ,w n ) that partitions C into regions of prescribed size is unique up 
to addition of the same constant to all W{. Consequently, the parts in the 
weighted Voronoi partition outside of C are also uniquely determined. 

The proof uses the following lemma. 

Lemma 2. Assume that C is path- connected, in addition to being an open 
bounded set. Consider a partition of the point set S in two nonempty sets 
T and S \ T. Suppose, that, for a given weight vector w, all regions C R 
VR w (p,S) have positive measure. 

If we increase the weight of every point in T by the same constant 5 > 0, 
then /i(C PI VR w (t, S)) increases, for some t G T. 

Proof. In this section, we will omit the reference to the set S from the Voronoi 
regions and simply write VR w (p) instead of VR w (p, S), since the point set S 
is fixed. 

It is clear that the regions VR u ,(t) for t G T can only grow, in the sense of 
gaining new points, but we have to show that this growth results in a strict 
increase of /x(C fl VR U) (t)) for some t G T. 

Consider the continuous function 

dr(z) := min{ d p (z) — w p | p G T } 
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and the function ds\r(z), which is defined analogously, and define 

VR W (T) := { z G R d | d T (z) < d s \ T (z) } 
VR W (S \ T) := { z G R d | d T (z) > d S \ T (z) } 

These two sets partition IR d , apart from some "generalized bisector" where 
equality holds. 

Roughly, VR W (T) is obtained by merging all Voronoi regions VR w (f) for 
t G T into one set, together with some bisecting boundaries between them. 
Since these boundaries have zero measure inside C, we have 

/i(CnVB 1fl (T)) = ^/i(CnVR tt (*)) (7) 

Let w = (wi, . . . , w n ) and w' = (w[, . . . , w' n ), where w[ := Wi + 5 if s, G T 
and w[ := otherwise. It suffices to show that in going from w to w', 
yu(C n VR„(T)) increases, since then, by ([7]), one of the constituent Voronoi 
regions ^(CnVR^t)) must increase. As mentioned, the set C nVR w (T) can 
only grow, but we have to show that this growth results in a strict increase 
of fl. 

If CnVR w '(S\T) = 0, we are done because /x(CflVR^/(T)) has increased 
to /i(C). Otherwise, consider a path in C connecting some point pi G VR W (T) 
with some point P2 G VR^/(5'\T). We have, by definition, dr(pi) < c?5\t(pi) 
and g?t(P2) — 5 > ds\r(P2)- Hence, by the intermediate value theorem, we 
can find a point p on the path with d T (p) —5/2 = ds\r(p)- This point is 
an interior point of C PI VR l „(S' \ T) and C fl VR U /(T), and hence there is a 
neighborhood of p which adds a positive measure to Cfl VR W (T) when going 
from w to w'. □ 

To finish the proof of uniqueness in Theorem [TJ suppose that two weight 
vectors w, w' do the job of producing the desired measures Aj for the regions, 
but w' — w (c, c, . . . , c) for any constant c. We may assume, by adding a 
constant to all entries, that w < w', and Wi = w\ for some i. We will now 
gradually change w into w' by modifying its entries. Let T := { s, G 5 | 
w ■ — = maxj(u^ — lUj) }. By assumption, T is neither empty nor equal 
to S. We increase the weights of the points in T, leaving the remaining values 
fixed. The amount 5 is chosen in such a way that w\ — Wi for i G T does not 
decrease below the remaining differences w'j — Wj for j G S\T. By Lemma |2j 
the measure of the region assigned to some point to £ T strictly increases in 
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this process. When the limiting value S is reached, the set T of points where 
w[ — Wi achieves its maximum is enlarged, and the whole process is repeated 
until w = w'. During this process, to will always be among the points whose 
weight is increased, and therefore, the measure of its region can never shrink 
back to the original value — a contradiction. □ 

6. Conclusion and future work 

We have given a geometric proof for the fact that additively weighted 
Voronoi diagrams can optimally solve some cases of the Monge-Kantorovich 
transportation problem, where one measure has finite support. Surprisingly, 
the existence of an optimal solution — the main mathematical challenge in 
the general case — is an easy consequence of our proof. In remains to be seen 
to which extent our assumptions on the distance functions d Pi can be further 
generalized. 
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